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^O , A subordinate Brownian motion is a Levy process which can be obtained by replacing the time of the 

Brownian motion by an independent subordinator. The infinitesimal generator of a subordinate Brownian 
motion is — 0( — A), where <j> is the Laplace exponent of the subordinator. In this paper, we consider a large 
class of subordinate Brownian motions without diffusion component and with <j> comparable to a regularly 
varying function at infinity. This class of processes includes symmetric stable processes, relativistic 
stable processes, sums of independent symmetric stable processes, sums of independent relativistic stable 
processes, and much more. We give sharp two-sided estimates on the Green functions of these subordinate 
Brownian motions in any bounded re-fat open set D. When D is a bounded C 1 ' 1 open set, we establish 
an explicit form of the estimates in terms of the distance to the boundary. As a consequence of such 
sharp Green function estimates, we obtain a boundary Harnack principle in C ' open sets with explicit 
rate of decay. 

AMS 2010 Mathematics Subject Classification: Primary 60J45, Secondary 60J75, 60G51. 
Keywords and phrases: Green function, Poisson kernel, subordinate Brownian motion, K-fat open set, 
C ' open set, symmetric stable process, relativistic stable process, harmonic functions, Harnack inequality, 
I/"") ' boundary Harnack principle, fluctuation theory, regularly varying function. 

in 

^ ! 1 Introduction 

The investigation of fine potential-theoretic properties of discontinuous Markov processes in the Euclidean 
space began in the late 1990's with the study of symmetric stable processes. One of the first results obtained 
in this area was sharp Green function estimates of symmetric a-stable processes in bounded C • domains 
inR d ,0<a<2,d>2. Recall that if A is a symmetric Markov process in R d and D is an open subset of 
R d , then the Green function Gd{x, y) of X in D (if it exists) is the density of the mean occupation time for 
H \ X before exiting D, that is, the density of the measure 
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U^E X lu(X t )dt, UcD, 

Jo 



where m is the first time the process X exits D. Analytically speaking, if L is the infinitesimal generator 
of X and C\d is the restriction of C to D with zero exterior condition, then Gd(',v) is the solution of 

{C\ D )U = ~Sy. 

A process X = (X t : t > 0) is called a (rotationally) symmetric a-stable (Levy) process, < a < 2, if 
it is a Levy process whose characteristic exponent $, defined by E[exp{i6* • X t }} — exp{— £$(#)}, is given 



'Supported by Basic Science Research Program through the National Research Foundation of Korea(NRF) grant funded by 
the Korea government(MEST)(2010-0001984). 

tSupported in part by the MZOS grant 037-0372790-2801. 



by Q(0) = \9\ a . The infinitesimal generator of a symmetric a-stable process is — (— A) a / 2 . The paths of 
the symmetric a-stable process X are purely discontinuous, as opposed to the case a = 2 corresponding 
to Brownian motion which has continuous paths. It was independently shown in [TT] and [35] that if D is 
a bounded C 1 ' 1 domain in R d , Gn(x,y) the Green function of the symmetric a-stable process in D, and 
5d{x) the distance between the point x and the complement D c of D, then there exists a constant c > 1 
(depending on D and a) such that 

--^..(MfQMiigy i <Ga( , < e f 1A fe(^fa))° / ') i i , (L1) 



for all x,y € D. Here and in the sequel, for a, 6 €E R, a A b := min{a, 6} and a V b := max{a, b}. The same 
form of the estimates in the case a = 2 (and d > 3) were obtained much earlier in [33] and |35j for the 
Brownian motion case. 

The proofs of (jl.ip for symmetric a-stable processes relied heavily on the explicit formulae for the Green 
functions and the Poisson kernels of the ball. Moving away from stable processes, such formulae were not 
available and new methods had to be developed. [25] studied the relativistic a-stable process (with relativistic 
mass m > 0) whose characteristic exponent is given by $(0) = (\8\ 2 +m 2 / a ) a / 2 — m and infinitesimal generator 
is given by m — (—A + to 2 / q ) q / 2 , and showed that the Green function of this process in any bounded C 1,1 
domain D satisfies the same sharp estimates (|1.1[) . Soon after, [12], using a perturbation method, established 
a general result which includes the main result of 25 as a special case. For different generalizations of the 
main result of [25] , see the recent papers [HI [19] . 

Quite recently, [8 j studied the Levy process which is the sum of independent symmetric /3-stable and 
a-stable processes, < j3 < a < 2. The characteristic exponent of this process is given by $(0) = \6\ a + \6\P 
and the infinitesimal generator by —(—A)"/ 2 — (— A)^/ 2 . Sharp two-sided estimates on the heat kernel of 
this process in C 1,1 open sets were established in [8]. As a by-product of the heat kernel estimates, sharp 
Green function estimates of the process in any bounded C ' open set were obtained in [8]. Estimates have 
the form (jTTTJ). In contrast with the relativistic stable processes, these Green function estimates cannot be 
obtained using the methods of [TJ] [TJ] [TJ] [35]. The case a = 2 (i.e., one of the processes is a Brownian 
motion) was covered in [10] with analogous estimates. 

The common feature of these Green function estimates is that both the distance between the points, 
\x — y\, and distances to D c , 8d(x),5d(ij), appear as arguments of power functions. However, it follows from 
[S] Chapter 5] that the asymptotic behavior of the free Green function G(x, y) of many transient symmetric 
Levy processes is of the form 

G ( x >y) ~ i \d- a oi\ pit as \ x ~ y\^° 

\x-y\ d a £{\x-y\ A ) 

where a € (0,2) and £ is a nontrivial slowly varying function at infinity. (See, also Theorem 12.91 below.) 
Therefore, Green function estimates of the form (jl.l[) cannot be true for these general symmetric Levy 
processes. The purpose of this paper is to establish sharp two-sided Green function estimates for these more 
general Levy processes in open sets of M. d . In our estimates, 5i){x),5D{y) and \x — y\ appear as arguments 
of regularly varying functions, not necessarily power functions. In order to explain our setting and results, 
let us first note that stable processes, relativistic stable processes and sums of independent stable processes 
can be obtained as subordinate Brownian motions. Indeed, let W = (W t = (W^,...,Wf) : t > 0) be a 
d-dimcnsional Brownian motion, and let S — (St : t > 0) be an independent subordinator. Recall that 
a subordinator is an increasing Levy process on [0,oo), which can be characterized through its Laplace 
exponent 4>: E[exp{-AS t }] = exp{-t^(A)}, A > 0. The process X = (X t : t > 0) defined by X t := Ws t is 
called a subordinate Brownian motion. The infinitesimal generator of X is — (f>(— A). By choosing the Laplace 
exponent 4>(X) as A"/ 2 , (A + m 2 / a ) a / 2 — m and A"/ 2 + A^/ 2 respectively, the resulting subordinate Brownian 



motion turns out to be a symmetric a-stable process, a relativistic stable process and an independent sum 
of j3 and a-stable processes respectively. The Laplace exponent of a subordinator is a Bernstein function 
and hence has the representation 

0(A) =b\+ I (1 ~ e- xt ) fi(dt) , 

J(0,oo) 

where b > and /i is a measure (called the Levy measure of 0) such that f,„ Jl A t) jj,{dt) < oo. If 
the measure /i has a completely monotone density, the Laplace exponent is called a complete Bernstein 
function. The common feature of the Laplace exponents 0(A) = A Q / 2 , cj)(X) = (A + m 2 / Q ) a / 2 — m and 
0(A) = A* 3 / 2 + A Q / 2 is that all three of them are complete Bernstein functions whose behavior at infinity is 
given by lim^oo 0(A)/A a / 2 = 1. We will see that those two properties (the latter slightly weakened) are the 
determining factors for the Green function estimates (jl.ljl . 

Recall that an open set D in R d (d > 2) is said to be a C 1,1 open set if there exist a localization radius 
R > and a constant A > such that for every z £ 3D, there exist a C^-function ip = ip z : E d_1 — > R 
satisfying ip(Q) = 0, V^(0) = (0, ...,0), ||Vt/>||oo < A, |V^(x) - V-0(z)| < A|x - z\, and an orthonormal 
coordinate system CS Z : y = (j/i, • • • , j/d-i, j/d) := (y> 2/d) with origin at z such that 

B(z, R)nD = {y= (y, y d ) £ B(0, R) in CS Z : y d > ^(y)}. 

The pair (R, A) is called the characteristics of the C 11 open set D. We remark that in some literature, the 
C 1 ' 1 open set defined above is called a uniform C 1 ' 1 open set since (R, A) is universal for all z £ dD. By a 
C 1 ' 1 open set in K we mean an open set which can be written as the union of disjoint intervals so that the 
minimum of the lengths of all these intervals is positive and the minimum of the distances between these 
intervals is positive. Note that a bounded C 1 ' 1 open set can be disconnected. 

The main result of this paper is the following sharp Green function estimates. In the statement and 
throughout the paper we use notation /(£) x g(t) as t — > oo (resp. t — > 0+) if the quotient f(t)/g(t) stays 
bounded between two positive constants as t — > oo (resp. t —¥ 0+). 

Theorem 1.1 Suppose that X = (X t : t > 0) is a Levy process whose characteristic exponent is given by 
$(6») = <j){\6\ 2 ), 9 £ R d , where 

<j> : (0, oo) — > [0, oo) is a complete Bernstein function such that 0(A) x A Q ' £(A) , A — > oo , (1-2) 

a £ (0, 2 A d) and £ : (0, oo) — > (0, oo) is a measurable, locally bounded function which is slowly varying at 
infinity. When d < 2, we assume an additional assumption, see (|2.14|) . Then for every bounded C ' open 
set D in M. d with characteristics (R,A), there exists C\ = C\{diam{D), R, A,a,£,d) > 1 such that the Green 
function Gn(x,y) of X in D satisfies the following estimates: 

c -i ( 1 A (6 D (x)6 D (y)r/H(\x-y\- 2 ) \ 1 

1 V v^((M*))- 2 K((M2/))- 2 ) l* - v\ a ) e (\ x ~ y\^\ x - y\ d ~ a 

< OntovXCL (lA ^121^ ) Tn * - d . (1.3) 

Since the subordinate Brownian motion X is completely determined by the Laplace exponent <fi, one 
would expect that the above estimates can be expressed in terms of the function only. And indeed, an 
alternative form of (jl.31) reads as follows: There exists c > 1 such that 

c -r( lA *(l* -Vl- 9 ) > \ 



V<t>(8n(x)- 2 )(t>(6 D (x)- 2 ) J \x ~ y\ d 0(\x - y\-i) 



<t>(\x-y\- 2 ) 
^/W^W^WJ^W^) J \x-y\ d <t>{\x-y\ 



< G D ( X ,y)<c[iA ll/r rrjL\, ,j ■- ■■„^_„ l -^ m 



(see (|1.6p below for yet another alternative form of these estimates). To the best of our knowledge, the above 
Green function estimates include all known Green function estimates of pure-jump transient subordinate 
Brownian motions in bounded C 1 ' 1 open set D in R d as special cases. The estimates above include a lot 
more processes: In the case d > 3, our estimates are valid for all subordinate Brownian motions satisfying 
(|1.2[) . For more concrete examples, see Example 12.161 

Let us give the main ingredients of the proof of Theorem 11.11 The groundwork has been laid down in 
the recent paper |20j where a similar class of subordinate Brownian motions was studied. One difference 
to the current setting was that in |20j the Laplace exponent was assumed to be precisely regularly varying 
at infinity and not just comparable to a regularly varying function. Another difference is that |20j contains 
some additional assumptions that we showed to be redundant. The results of [3U] are reproved in [3T] under 
conditions valid in this paper (with the redundant assumptions removed). When referring to those results 
we will quote both sources. The main result of [20[ 121) is the boundary Harnack principle for nonnegative 
harmonic functions of the subordinate Brownian motion X in bounded K-fat open sets. Based on the 
boundary Harnack principle and using the well-established methodology of 4, 16 , we will first obtain Green 
function estimates of the form (|1.5p (in the spirit of [4] [16]) in bounded K-fat open sets. 

Recall from [32] that an open set D in M. d is K-fat if there exists R\ > such that for each Q £ 3D and 
r £ (0, i?i), D n B(Q, r) contains a ball B(A r (Q), nr). The pair {R\, k) is called the characteristics of the 
K-fat open set D. All Lipschitz domains, non-tangentially accessible domains and John domains are K-fat 
(cf. [18l|32] and the references therein). In general, the boundary of a K-fat open set can be nonrectifiable. 

Theorem 1.2 Suppose that X = (X t : t > 0) is a Levy process satisfying the same conditions as in 
Theorem, \1.1\ and that D is a bounded K-fat open set with characteristics (Ri,k). Then there exists Ci = 
C2(diam(D), Ri, K,a,£,d) > 1 such that for every x, y £ D, 

n -i g(x)g(y) ^ n , \ <, n 9{x)g{y) , _ „, , n K , 

62 g(A)*\x - y\*-°t{\z - y\-*) ~ G °^ V) ~ ° 2 g(A^\x - y\*-»t{\x - y\^) ' i4e ^ a! 'W' (L5) 
where g and B(x, y) are defined in (J3.11I) and (|3.7p respectively. 

In the case < C\ < l(X) < C2 < 00 for large A, using the Harnack inequality and the boundary Harnack 
principle, the above form of Green function estimates has been established by several authors in special 
cases. See [TUl Theorem 1.1], [13 Theorem 2.4] and [T7] Theorem 1]. 

To obtain the interior estimates in Theorem [L2] (i.e. for points x, y away from the boundary), we use the 
asymptotic behavior of the Green function of X in M. d proved in [21] Theorem 3.2] (see also [20] Theorem 
3.1]). Using the interior estimates and following the method developed in [4] [16], we obtain the full estimates 
in a bounded K-fat open set using the boundary Harnack principle from [20[ 121] . 

Even though a lot of complications occur due to the appearance of the slowly varying function £, the 
proof of Theorem 11.21 is still routine. But the precise estimates (|1.3|) in bounded C 11 open sets are very 
delicate. One of the ingredients comes from the fluctuation theory of one-dimensional Levy processes. Let 
Z = (Z t : t > 0) be the one-dimensional subordinate Brownian motion defined by Z t :— Wg , and let V be 
the renewal function of the ladder height process of Z. The function V is harmonic for the process Z killed 
upon exiting (0, 00), and the function w(x) := U(xd)l/ Xd>0 }, x G K d , is harmonic for the process X killed 
upon exiting the half space Ri := {x = (x\, . . . , Xd-i, Xd) £ R d : Xd > 0} (Theorem 14. ip . Therefore, w gives 
the correct rate of decay of harmonic functions near the boundary of Ri. This shows the importance of the 
fluctuation theory (of one-dimensional Levy processes) in our approach. 



The second ingredient is the "test function" method applied to the operator A defined by 
Af (x) = lira [ (f(y)-f(x))J(y-x)dy, 

e ^" J{yeWl d :\y-x\>e} 

with the domain consisting of functions / for which the limit exists and is finite. Here J denotes the density of 
the Levy measure of X. On the space of smooth functions with compact support, this operator coincides with 
the infinitesimal generator of X. We emphasize that, compared to the test function methods of [51 1§1 [T5]. 
there are several differences in our approach. In [5J |H1 115) . appropriate subharmonic and superharmonic 
functions of X (or the truncated version of X) are chosen as test functions, first in the case of half spaces 
and then for C ' open sets, and the values of the generator acting on these test functions are computed 
in detail. Then suitable combinations of the test functions are used to find the correct exit distribution 
estimates. In [5J [9j [T5], the test functions are power functions of the form x — > (xd) p and the densities of the 
Levy measures of the processes have closed forms. However, the density J of the Levy measure of our process 
does not have a simple form. We do not even know the asymptotic behavior of J near infinity in general. 
Furthermore, in general, power functions of the form x — > (xd) p are neither subharmonic nor superharmonic 
functions for our processes, and it is not clear what are the appropriate choices for the test functions. 

Due to the above differences and difficulties, obtaining the correct boundary decay rate of the Green 
function in C ' open set D requires new ideas and approaches. In this paper, we will use the function w 
which is smooth and harmonic on the half space, as our only test function. Using this and the characterization 
of harmonic functions recently established in [7], we show that Aw = on the half space (Theorem 14. 3|) . 
With this, we prove the following fact in Lemma 14.41 which is the key to the proof of Theorem 11.11 If D 
is a C 1 ' 1 open set with characteristics (i?,A), Q £ dD and h(y) = V(5D{y))lDnB(Q,R), then Ah(y) is well 
defined and bounded for y £ D close enough to the boundary point Q. Using this lemma, we give certain 
exit distribution estimates in Lemma 14. 5[ which provide the correct rate of decay of Green functions near 
the boundary of D. Unlike El [15], in Lemma F4.5I we do not construct subharmonic and superharmonic 
functions on C 1 ' 1 open set D. Instead we use Dynkin's formula on h to obtain the desired exit distribution 
estimates directly. In fact, our approach is simpler than the previous approaches and may be used for other 
types of jump processes. We hope our approach will shed new light on the understanding of the boundary 
behavior of nonnegative harmonic functions of general Markov processes. 

The estimates (|1 .3[) are best understood in terms of the renewal function V which provides the exact rate 
of decay of Gd near the boundary. Let G be the Green function of X in the whole space M. d . An equivalent 
form of (|1.3[) is given by 

c [ 1A V{\x-y\Y )G(x,y)<G D (x,y)<c{lA v{]x _ y ^ ) G(x,y). (1.6) 

By combining the sharp estimates of the Green function in a bounded C 1,1 open set with the boundary 
Harnack principle proved in |20U21] (see Theorem [2T5J below), we obtain a boundary Harnack principle with 
explicit decay rate. In the next theorem we give an extension to unbounded C ' open sets. Recall that, 
given Q £ dD, a function u : M. d — \ M. is said to vanish continuously on D c C\B(Q,r) if u — on D c P\B(Q, r) 
and u is continuous at every point of dD n B(Q, r). 

Theorem 1.3 Suppose that X = (X t : t > 0) is a Levy process satisfying the same conditions as in Theorem 
\1.1\ and that D is a (possibly unbounded) C 1 ' 1 open set in M. d with characteristics (R,X). Then there exists 
C*3 = Cs(R, A, a, £,d) > such that for r £ (0, (R A l)/4], Q £ dD and any nonnegative function u in E d 
that is harmonic in D n B(Q, r) with respect to X and vanishes continuously on D c n B(Q, r), we have 

" ir] < c 3 -— — ^S==== for every x,y€DnB(Q,r/2). (1.7) 



M*) Q/ V^((M2/))- 2 ) " M2/) Q/ V^((M*))- 



An alternative form of (|1.7[) reads as follows: There exists a constant c > such that 

" ! ' < c - "^ .. for every x,y e D D B(Q,r/2). (1.8) 



y(fo(x)) " V{8 D {y)) 

Note that unlike the usual form of the boundary Harnack principle where one considers the ratio of two 
harmonic functions, functions in the denominator of (jl.71) and (jl.8p are not harmonic. Instead, they provide 
the correct boundary decay of non- negative harmonic functions. Indeed, an equivalent form of Theorem 11.31 
says that there exists a constant c > 1 such that for any nonnegative function u in M. d that is harmonic in 
D n B(Q, r) with respect to X and vanishes continuously on D c n B(Q, r) it holds that 

c~ WiSoix)) < u(x) < cV(S D (x)) for every leflfl B(Q, r/2). 

This paper is organized as follows: In the next section we establish the setting and notation, prove 
several new results for complete Bernstein functions, and describe some of the known results from |20[ 121) . 
In Section 3 we prove the Green function estimates in bounded /c-fat open sets. Section 4 is devoted to the 
Green function estimates in bounded C 1 ' 1 open sets. 

We will use the following conventions in this paper. The values of the constants C\, C2, ■ ■ ■ , M, e\ and 
R, R\, i?2, • • • will remain the same throughout this paper, while c, Co, c\, C2, • • • and r, ro, r±, r<i, . . . stand 
for constants whose values are unimportant and which may change from one appearance to another. All 
constants are positive finite numbers. The labeling of the constants Co, Ci, C2, • • • starts anew in the statement 
of each result. The dependence of the constants on dimension d may not be mentioned explicitly. We will 
use ":=" to denote a definition, which is read as "is defined to be". Further, f(t) ~ g(t), t — > (/(£) ~ g(t), 
t — > 00, respectively) means lim t _s.o f(t)/g(t) — 1 (lim f _i. 0O f(t)/g(t) = 1, respectively). For any open set U, 
we denote by 5jj{x) the distance between x and the complement of U, i.e., Sjj(x) = dist(ic, U c ). We will use 
d to denote the cemetery point and for every function /, we extend its definition to d by setting f(d) = 0. 
For every function /, let /+ := / V 0. We will use dx to denote the Lebesgue measure in M d . For a Borel set 
A C R d , we also use \A\ to denote its Lebesgue measure and diam(A) to denote the diameter of the set A. 

2 Preliminaries 

In this section we collect and explain preliminary results necessary for further development in Sections 3 
and 4. Most of these results originate from [20) where they were proved under somewhat stronger conditions 
than in this paper. Their extensions to the current setting, in particular Theorems 12.91 12.111 and 12.151 are 
given with full proofs in |21j . Here we prove only results that have not appeared in |20| . Lemma 12.11 and 
Propositions 12.41 and 12.61 about complete Bernstein functions may be of independent interest. The difference 
between the assumptions in [5D] and this paper is discussed in Remark 12.21 

A C°° function <f> : (0, 00) — > [0, 00) is called a Bernstein function if (—l) n D n <fi < for every positive 
integer n. Every Bernstein function has a representation </>(A) = a + bX + J, Q Jl — e~ xt ) n{dt) where 
a, b > and /1 is a measure on (0, 00) satisfying J, Q >(1 A t) n(dt) < 00; a is called the killing coefficient, 
b the drift and /x the Levy measure of the Bernstein function. A Bernstein function </> is called a complete 
Bernstein function if the Levy measure \x has a completely monotone density n(t), i.e., {—l) n D n ^ > for 
every non-negative integer n. Here and below, by abuse of notation we will denote the Levy density by fi(t). 
For more on Bernstein and complete Bernstein functions we refer the readers to [27] . 

First, we show that <f> being a complete Bernstein function implies that its Levy density cannot decrease 
too fast in the following sense: 

Lemma 2.1 Suppose that <fi is a complete Bernstein function with Levy density fi. Then there exists C4 > 1 
such that /i(t) < Cifi(t + 1) for every t > 1. 

6 



Proof. Since /i is a completely monotone function, by Bernstein's theorem ( |27| Theorem 1.4]), there 
exists a measure m on [0, oo) such that /x(i) = J, Q , e~ tx m(dx). Choose r > such that J, Q , e _a: m(dx) > 
J, s e _a: m(dx). Then, for any t > 1, we have 

e~ t:,; TO(dx) > e - ( *- 1)r / e- x m(d^) 

[0,r] ^[0,r] 



> e-C^ 1 ^ / e - x m(dx) > / e~ tx m{dx). 

J (r,oo) J(r,oo) 



Therefore, for any t > 1, 



M(t + 1) > / e^ 1 ^ m(dx) > e- r / e"** m(cfe) >le~ r [ e' tx m(dx) = ^ e - r /i(t). 

J[0,r] J[0,r] 2 J[0,oo) 2 



□ 



Suppose that S = (St : t > 0) is a subordinator with Laplace exponent </), that is 

E[e~ ASt ] =e- l ^ X) , Vt,A>0. 

The Laplace exponent of a subordinator is always a Bernstein function. Let C/(A) := E L lig t £A} dt 
denote the potential measure of S. If is a complete Bernstein function with infinite Levy measure, then 
the potential measure U has a completely monotone density u(t) (see, e.g., (27| Remark 10.6 and Corollary 
10.7]). 

Recall that a function £ : (0, oo) — > (0, oo) is slowly varying at infinity if 

£(Xt) 
lim — 1 , for every A > . 

In the remainder of this paper we assume that is a complete Bernstein function and we will always 
impose the following 

Assumption (H): There exist a € (0,2) and a function I : (0,oo) — > (0,oo) which is measurable, locally 
bounded and slowly varying at infinity such that 

^(A)xA a/2 l(A), A^oo. (2.1) 

Remark 2.2 (a) The precise interpretation of (j2.ll) will be as follows: There exists a positive constant c > 1 
such that 

c_1 " a^aT - c **«**e[i,oo). 

The choice of the interval [l,oo) is, of course, arbitrary. Any interval [a,oo) would do, but with a different 
constant. This follows from the assumption that £ is locally bounded. Moreover, by choosing a > large 
enough, we could dispense with the local boundedness assumption. Indeed, by [H Lemma 1.3.2], every slowly 
varying function at infinity is locally bounded on [a, oo) for a large enough. Although the choice of interval 
[1, oo) is arbitrary it will have as a consequence the fact that all relations of the type /(£) x g(t) as t — > oo 
(respectively t — > 0+) following from (|2.ip will be interpreted as c~ l < f(t)/g(t) < c for t > 1 (respectively 
< t < 1) for an appropriate constant c. 

(b) The assumption (H) is an assumption about the behavior of </> at infinity. We make no assumption on 
</> near zero. As a consequence, we will be able to obtain information about the small scale behavior of the 
subordinate process, but almost nothing can be inferred about its large scale behavior. 

7 



(c) The main assumption in |20j was that is a complete Bernstein function such that 

</>(A) = \ a/2 £{\) , for all A > , (2.2) 

where a £ (0, 2) and £ is a slowly varying function at infinity. This assumption allows us to obtain exact 
asymptotic behavior of various functions. More precisely some of the results in }20| were of the form 
f(t) ~ g(t), while with the assumption (12.1[) we can obtain only the corresponding results in the weaker 
form f(t) x g(t). Proofs of these weaker results can be found in [21]. We note that our current assumptions 
are indeed strictly weaker than the ones in [20]: There exists a complete Bernstein function satisfying (|2.1[) 
which is not regularly varying at infinity, see [21] Example 2.8]. 

(d) We briefly comment on the other assumptions from [20j which are now removed. The assumption Al 
in [20] needed for transience in case d < 2 is replaced by (|2.14|) below. The assumptions A2 and A3 in [20] 
used for the technical lemma [6] Lemma 5.32] are redundant - see [2"T] Lemma 3.1]. Assumption A4 in [20] 
is always valid for complete Bernstein function as proved here in Lemma 12. II Finally, the assumption (2.5) 
in [20] Proposition 2.2] is no longer needed as [20[ Proposition 2.2] is now replaced by Proposition ^. 6l below. 

It follows from (|2.1j) that liniA-s-oo 4>(X)/X = and lim>_ ! . oc 4>(X) = oo, implying that (j> has no drift and 
its Levy measure is infinite. Therefore, the potential measure U of the corresponding subordinator S has a 
completely monotone density u. 

The behavior of u(t) and the density /i(t) of the Levy measure can be inferred from the following result. 

Proposition 2.3 (34, Theorem 7]) Suppose that ip is a completely monotone function given by 

V>(A)=/ e- xt f(t)dt, 
Jo 

where f is a nonnegative decreasing function. Then 

f(t)<(i-e- l y l t- l ^(r l ), t>0. 

If, furthermore, there exist 5 £ (0, 1) and a, to > such that 

ip(r\) <ar" 5 V(A), r>l,t>l/t , (2.3) 

then there exists C$ = C$(w, f,a,to,S) > such that 

/(*) > Csrvr 1 ), t < t . 



We first apply the above proposition to ip(X) — <fi(\) 1 — J Q e xt u(t) dt to obtain the behavior of u near 



zero: 



^xrVrVx^, t->0+. (2.4) 

Condition (|2.3[) follows from (|2.1[) by use of Potter's theorem (cf. [2] Theorem 1.5.6]). By applying (|2.4[) to 
the complete Bernstein function A n- X/cj>(X) ([27] Proposition 7.1]) one obtains the following behavior of 
/j,(t) near zero: 

Ktjxi-^r^r^r 1 ), t^o+. (2.5) 

We refer the reader to [2D Theorem 2.9, Theorem 2.10] for the detailed proofs of (j2~4")l and ([2~5]) . The 
corresponding precise asymptotics are given in [23 p. 1603] under the assumption (|2.2I) . 

A consequence of the asymptotic behavior (|2.5[) of /i(t) is that for any K > there exists c = c(iT) > 1 
such that 

fi(t) < cfi(2t), Vte(0,K). (2.6) 



The behavior of /j,(t) at infinity has already been determined in Lemma |2. II There exists a constant c > 1 
such that 

fi(t)<cn(t+l), VOL (2.7) 

This property of fj, was assumed in [20] as A4, but we have shown in Lemma l2.1l that it always holds true. 

We consider now one-dimensional subordinate Brownian motions. Let B = (B t : t > 0) be a Brownian 
motion in K., independent of S 7 with 



E 



je(B t -B ) 



-W 



V6»eK,i>0. 



The subordinate Brownian motion Z = (Zt : t > 0) in R defined by Zt = _Bs t is a symmetric Levy process 
with the characteristic exponent $(6>) = <p(0 2 ) for all 9 £ R. 

Let Z t := sup{0 V Z 8 : < s < t} be the supremum process of Z and let L = (L t : £ > 0) be a local time 
of Z — Z at 0. i is also called a local time of the process Z reflected at the supremum. The right continuous 
inverse L~[ of L is a subordinator and is called the ladder time process of Z . The process Z L -i is also a 
subordinator and is called the ladder height process of Z. (For the basic properties of the ladder time and 
ladder height processes, we refer our readers to [TJ Chapter 6].) 

Let x be the Laplace exponent of the ladder height process of Z. It follows from [T3] Corollary 9.7] that 



m ( i r io s' 



*SPU)--P(i 



log(<KA 2 fl 2 )) 

1 + 2 



rf(9 



VAX). 



(2.8) 



Using (|2.8j) . it was proved in [20j Proposition 2.1] that, if is a special Bernstein function, so is Xj i- e -j 
A h^ A/x(A) is also a Bernstein function. The next result tells us that such relation is also true for complete 
Bernstein functions. For the proof we will need the following fact, see (27j Theorem 6.10]: If <f> is a complete 
Bernstein function, then there exist a real number 7 and a [0, l]-valued function r\ on (0, 00) such that 



log 0(A) =7 + 



1 + 1 2 X + t 



r)(t)dt. 



(2.9) 



Proposition 2.4 Suppose <f>, the Laplace exponent of the subordinator S , is a complete Bernstein function. 
Then the Laplace exponent \ °f the ladder height process of the subordinate Brownian motion Z t = B$ t is 
also a complete Bernstein function. 



OO /-0O 



Proof. By (J2T8J and (|2l)]) , we have 

io gX (A)=j+- r 1 

2 7T Jo Jo 
By using < n(t) < 1, we have 



Since 



n{t) 



t 



1 



1 + 1 2 x 2 e 2 + 1 



■ de = - 



t 



1 



< 



< 



1 + t 2 X 2 6 2 + t 



1 



1 + t 2 1 
1 



1 



n(t)dt- 



cW 



1 



X 2 9 2 t 



x 2 e 2 + t x 2 e 2 + t 
x 2 t 



1 + t 2 \ x 2 e 2 + 1 x 2 e 2 + 1 



■ d* = i^ 



»o± + 1 t A 7o 7 2 + 1 



c?7 : 



2A v / i' 



we can use Fubini's theorem to get 



logx(A) 



2 

2 " 

7 
2 ' 

7i 



t 



1 



2(1 + t 2 ) 2Vt(X + Vt) 
t 1 



r)(t)dt 



(2.10) 



2(1 + t 2 ) 2(1 + 1) 
s 1 



»7(t)d!< 



\2(l + t) 2V~t{\ + V~t) 



r)(t)dt 



r)(s )ds, 



□ 



/o \l + s 2 A + s, 
Applying 27, Theorem 6.10] we get that \ is a complete Bernstein function. 

Remark 2.5 The above result has been independently proved in [23l Lemma 4]. 

The next result relates the behavior of \ with that of <j>. It will be used to obtain the asymptotic behavior 
of x at infinity. 

Proposition 2.6 Suppose that (f>, the Laplace exponent of the subordinator S, is a complete Bernstein 
function. Then the Laplace exponent \ of the ladder height process of Z satisfies 



e-"/ 2 v ^(A2)<x(A)<e" /2 v / 0(A2), for all X > . 

Proof. By the representations (|2.9|) and (|2.10|) . we get that for all A > 



(2.11) 



logx(A)--log<MA 2 ) 



< 



1 + t 2 Vt(A + Vt) 

\(Vt + x) 



t 



1 



■dt = 



1 + 1 2 X 2 + t 
X 



rj(t) dt 



dt = 



2 J {X 2 +t)^t(X + V~t) 2 7 (X 2 + t)V~t ' 2 

This implies that -vr/2 < logx(A) - ilog0(A 2 ) < tt/2 for every A > 0, which is (|2~TT|) 



□ 



Remark 2.7 We note that for the last two propositions we only need to assume that <j> is a complete 
Bernstein function; the assumption (H) is not used. 

Let V denote the potential measure of the ladder height process of Z. We will also use V to denote the 
corresponding renewal function, V(t) := V((0,t)). It follows from (|2.1|) and (|2.11|) that liniA^oo x(A)/A = 
and lim>,->oo x(A) = oo. Therefore, the ladder height process of Z has no drift and has infinite Levy measure. 
This suffices to conclude that the potential measure V has a density denoted by v, and the renewal function 
can be written as V(t) — L v(s) ds. Since x is a complete Bernstein function, v is completely monotone. 
We record these facts as 

Corollary 2.8 Suppose <f>, the Laplace exponent of the subordinator S , is a complete Bernstein function 
satisfying the assumption (H). Then the potential measure of the ladder height process of the subordinate 
Brownian motion Z t = B$ t has a completely monotone density v. In particular, v and the renewal function 
V are C°° functions. 

The smoothness of the renewal function V of the ladder height process Z will be used later in this paper. 

Similarly to the case of the density u of the potential measure U of the subordinator S in (|2.4[) . by using 
Proposition 12.61 we can obtain the asymptotic behavior of the renewal function V and its density v of the 
ladder height process of Z: 

v(t) x 0r 2 r 1/2 x {e{ p ))1/2 , ^o + , (2.i2) 
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v(t) x r^ft- 2 )- 1 / 2 x {l ^ 1/a , t^0+, (2.13) 

see [211 Proposition 3.9]. The corresponding precise asymptotics are given in 20, Proposition 2.7] under the 
assumption (J2.2I) . 

We next consider multidimensional subordinate Brownian motions. Let W = (Wt = (W/, . . . , W d ) : t > 
0) be a Brownian motion in M. d with 



E 



e i$-(W t -Wo) 



e~ m , V(9eIR d ,£>0, 



and let S be a subordinates independent of W with Laplace exponent 4>. In the remainder of this paper we 
will use X — (X t : t > 0) to denote the subordinate Brownian motion defined by X t = Ws t ■ The process 
X is a (rotationally) symmetric Levy process with the characteristic exponent given by &(9) = cf>(\8\ 2 ). It is 
easy to check that when d > 3 the process X is transient. This follows from the criterion of Chung- Fuchs 
type (e.g., [551 P- 252]) which for the subordinate Brownian motion X translates to the following: X is 
transient if and only if 

Liw dX<+cc - 

Since transience is a global property of the process, it cannot be inferred from the behavior of <p at infinity. 
For example, 4>(X) = log(l + A) + A"/ 2 , a € (0, 2), is a complete Bernstein function satisfying (j2.ll) . but the 
corresponding subordinate Brownian motion is recurrent in dimensions 1 and 2. To ensure transience, we 
will assume that in the case d < 2, there exists 7 £ [0, d/2) such that 

liminf^>0. (2.14) 

An immediate consequence of this assumption and [21] Corollary 2.6] is that the potential density u of S 
satisfies u(t) < ct 1 ^ 1 for all t > 1, where c > is some positive constant (cf. assumption Al from [20]). 

Transience of the process X ensures that the Green function G(x, y), x, y € M. d , is well defined. By spatial 
homogeneity we may write G(x,y) — G(x — y), where the function G is radial and given by the following 
formula, 

G(x) = / (iTTt)- d / 2 e-^ 2 /^u{t)dt 7 x G R d . (2.15) 

Jo 

Since u is decreasing, we see that G is radially decreasing and continuous in M d \ {0}. 
The Levy measure of the process X has a density J, called the Levy density, given by 

/>oo 

J{x)= {ATTt)- d l 2 e-W 2 / { ^ii{t)dt, a;eK d . 

Jo 

Thus J(x) = id^D with 

/•oo 

j(r):= (47rf)- rf/2 e- r2/(4t V(*) dt, r>0. (2.16) 

Jo 

Note that the function r 1— > j(r) is continuous and decreasing on (0, 00). We will sometimes use the notation 

J(x,y) for J(x-y). 

We discuss now the behavior of G and j near the origin. Under the assumption (|2.2[) . the precise 

asymptotic behavior was obtained in |201 Theorem 3.1, Theorem 3.2] by using precise asymptotic behavior 

of the potential density u, and, respectively, Levy density /i. These two results were proved by use of [5J 

Lemma 5.32], which required additional assumptions which were stated as A2 and A3 in 20 . It turned out 

that by using Potter's theorem (2, Theorem 1.5.6]) one can circumvent those assumption and still obtain 
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the conclusion of the lemma. The details are provided in [^TJ Lemma 3.1]. The lemma combined with (|2.4[) 
(resp. (|2.5jl ) and representation (|2.15[) (resp. (|2.16|l ) gives the following asymptotic behavior of G (resp. J 
under the assumption (I2.1[) : 

Theorem 2.9 ( |21[ Theorem 3.2]) Suppose that the Laplace exponent <f> is a complete Bernstein function 
satisfying the assumption (H) and that a € (0, 2 A d). In the case d < 2, we further assume (12.141) . Then 

G( - X ^~ \x\ d 0{\x\~ 2 ) ~ |a;| d -<^(|2;|- 2 )' '*' ~* °' 

Remark 2.10 Since a is always assumed to be in (0, 2), the assumption a € (0, 2 A d) in the theorem above 
makes a difference only in the case d = 1. 

Theorem 2.11 ( |21[ Theorem 3.4]) Suppose that the Laplace exponent <f> is a complete Bernstein function 
satisfying the assumption (H). Then 

Using (|2.6|) and (|2.7|) . and repeating the proof of [24j Lemma 4.2] we get that 

(a) For any K > 0, there exists c = c(if ) > 1 such that 

j(r)<cj(2r), Vr€(0,JO. (2.17) 

(b) There exists c > 1 such that 

j(r) < cj(r+l), Vr>l. (2.18) 

Therefore by [30l Theorem 2.2 and Section 3.1] (see also [6j Theorem 5.66], [211 Theorem 4.7, Corollary 4.8] 
and [24]) the Harnack inequality is valid for the process X. Before we state the Harnack inequality, we recall 
the definition of harmonic functions. 

For any open set D, we use td to denote the first exit time of D, i.e., td = mf{£ > : X t ^ D}. 

Definition 2.12 Let D be an open subset ofM. . A function u defined on R is said to be 

(1) harmonic in D with respect to X if E^ [|u(X TB )|] < oo and u{x) = E^ [u(X TB )], i£B, for every open 

set B whose closure is a compact subset of D; 

(2) regular harmonic in D with respect to X if it is harmonic in D with respect to X and for each x £ D, 

u(x)=E x [u{X TD )}. 

Theorem 2.13 (Harnack inequality) There exists C§ > such that, for any r € (0,1/4), xq £ R d , and 
any function u which is nonnegative on R and harmonic with respect to X in B(xo, 16r), we have 

sup u(y) < C 6 inf u{y). 

y£B(x ,r/2) y£B(x a ,r/2) 

^From now we will always assume that is a complete Bernstein function satisfying the assumption (H) 
for a £ (0, 2 A d) and the additional (J2.14J) in the case d < 2. We will no longer explicitly mention these 
assumptions 

For any open set D in R d , we will use Go(x,y) to denote the Green function of X in D. Using the 
continuity and the radial decreasing property of G, we can easily check that Go is continuous in (D x D) \ 
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{(x, x) : x £ D}. We will frequently use the well-known fact that Go (',2/) is harmonic in D\{y}, and regular 
harmonic in D \ B(y, e) for every e > 0. 

Using the Levy system for X, we know that for every bounded open subset D and every / > and 



E, [f(X TD ): X TD _ ± X TD ] = [_ [ G D (x lZ )J(z - y)dzf(y)dy. 



(2.19) 



<D JD 

Now we prove the following version of the Harnack inequality for X. 



Theorem 2.14 Let L > 0. There exists a positive constant Cj — Cj(L) > such that the following is true: 
If x\,xi £ K d and r £ (0, 1) are such that \xi — X2\ < Lr, then for every nonnegative function u which is 
harmonic with respect to X in B(x\,r) U B{x 2l r), we have 

Cj u(x 2 ) < u(x\) < Ctu{x 2 ). 

Proof. By [2TJ Proposition 4.10] (see also [501 Proposition 3.8]), for every r £ (0, 1), every lef 1 and every 
y £ B[x, §) it holds that 



K B (x^)(x,y) 



B(x,%) 



Gb{ x ,%){x,z)J{z - y)dz > cxj(\x-y\) 



£(r- 



(2.20) 



with some positive constant c\ > 0. 

Now let r £ (0, 1), Xi,x 2 £ K d be such that \xi — x 2 \ < Lr and let u be a nonnegative function which is 
harmonic in B(xi,r) U B(x2,r) with respect to X. 

If \xi — x 2 \ < ir, then since r < 1, the theorem is true from Theorem l2.13l Thus we only need to consider 
the case when |r < \xi — x 2 \ < Lr with L > |. 

Let u> S B(xi, §). Because |x 2 — Ztf| < |xi — x 2 \ + \w — Xi\ < (L+ |)r < 2Lr, by the monotonicity of j 
and f2~20| 



^B(i J ,?)(!t2i«') > cij(2Lr)- 



(2.21) 



7(r- 2 ) ' 
For any y £ B(x\, |), u is regular harmonic in B(y, ^) U S(xi, -y). Since \y — Xi\ < |, by Theorem 12.131 



u(y)>c 2 u(x 1 ), !/£%,-), 
for some constant c 2 > 0. Therefore, by (j2~il)]) and ([2T21"j) - (j2~2"2"]l . 



(2.22) 



u(x 2 ) = E X2 \u{X TB{x2i) ) 



> E,r 2 



> C 2 U 



(xi)P X2 (X T £%,r) =c 2 u(xi) / if B(cC2 r)(a;2,w)dw 

v ' s 8 ' Jb{x 1 ,d 



> cau(xi) B(xi,-) j(2Lr) 



c i u(x 1 )j{2Lr)- 



r a+d 



•£(r-*) ~*~ v 1,jy '£{r- 2 )' 

Thus, by Theorem 12.111 and the inequality above, there exists a constant C5 > such that for all r £ 
(0, 1), u(x 2 ) > c^u(xi) 3 .,J' . The right-hand side is by (|2.1T[) greater than C5C 6 ° S ' ° s u{x\) where c 6 = 
C z f 1 (4L) £ (0, 1). We have thus proved the right-hand side inequality in the conclusion of the theorem. The 
inequality on the left-hand side can be proved similarly. □ 



In [5D] we have established the boundary Harnack principle under the assumption (|2.2p (and the additional 
transience assumption in case d < 2) for K-fat open sets. Even though we only explicitly stated the results for 
d > 2 in [20) . the results and arguments there are in fact valid for d = 1 also. Under the current assumptions, 
the same result is reproved in [21) . 
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Theorem 2.15 ([501 Theorem 4.8], [HJ Theorem 4.22]) Suppose that I? is a bounded K-fat open set with 
the characteristics (Ri,n). There exists a constant C$ — Cg(d,a,l,Ri,K) > 1 such that if r < Ri A | and 
Q G dD, then for any nonnegative functions u, v in M. d which are regular harmonic in D n B(Q, 2r) with 
respect to X and vanish in D c n -B(Q, 2r), we have 



,_i «(A-(Q)) u(g) u(A-(Q)) t , nnBfo r ( 

8 ^4^Q)) "^) - C8 <I^Q))' l6DnB(ft 2' 



Before concluding this section, we give some examples satisfying our assumptions. 
Example 2.16 Suppose that < (3 < a < 2, 0<7<2 — a and define 

0i(A) = A a / 2 , 2 (A) = (A + l) Q / 2 -l, 03(A) = X a/2 + X^ 2 , 

4 (A) = A Q / 2 (log(l + A))^/ 2 and ^(A) = A Q / 2 (log(l + A))-- 9 / 2 . 
Then 0j, i = 1, . . . , 5, arc complete Bernstein functions which can be written as 

4 (A) = A Q/2 4(A), * = 1,...,5, 

with 

4(A) = 1, £ 2 (A) = ((A + l) a / 2 -l)A- Q / 2 , l s (\) = l + \V > - a V 2 , 

4(A) = (log(l + XW' 2 and 4(A) = (log(l + X))- fi ' 2 . 

As already mentioned in the introduction, the subordinate Brownian motion corresponding to 0i is a sym- 
metric a-stable process, the subordinate Brownian motion corresponding to 02 is a relativistic a-stable 
process and the subordinate Brownian motion corresponding to 03 is the sum of a symmetric a-stable pro- 
cess and an independent symmetric /3-stable process. The subordinate Brownian motions corresponding to 
04 and 05 were discussed in [6]. 

In the case d > 3, the only condition on the complete Bernstein function is (|1.2|) . so we can use 
Proposition 7.1, Corollary 7.9, Propositions 7.10-7.11, Corollary 7.12 of [27] to come up with infinitely many 
examples of such functions, e.g.: (i) A Q / 2 (log(l + log(l + A^'/ 2 )' 5 / 2 ))' 3 / 2 , aj.ie (0,2),/3 G (0,2 - a]; (ii) 
A Q / 2 (log(l + log(l + A^/ 2 ) d "/2))-/3/2 : a , 7 , 5 € (0, 2), /3 e (0, a]. 

All of the listed example satisfy the stronger condition (|2.2|) . As already mentioned, [2TJ Example 2.8] 
provides an example of a complete Bernstein function which satisfies (|2.1|) . but not (|2.2j) . 



3 Green function estimates on bounded K-fat open sets 

In this section we will establish sharp two-sided Green function estimates for X in any bounded K-fat open 
subset D of R d , d > 1. Our standing assumption is that is a complete Bernstein function satisfying the 
assumption (H) for a e (0, 2 A d) and the additional assumption (|2.14|) when d < 2. 

Lemma 3.1 There exist R2 — R^il) > 0, L\ > 2 and Cg > smc/j t/iai 

G(s) - G(L lX ) > C 9 i^id-^^^i-2) for every \x\ < R 2 . 
Proof. By Theorem 12 .91 there exists a constant c > 1 such that for all x € M d \{0} with \x\ < 1 it holds that 

1 —-< G(x) < 



\x\ d -»£(\x\- 2 ) ~ v ' ~ \x\ d -«l(\x\-' 1 ) ' 
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Choose Li — (4c 2 ) d -° V 2 so that c 2 /Lf a < 1/4. Since £ is slowly varying at infinity, there exists r\ < 1 

such that 

£(|a;|-^ 



t(\L lX \-i) 



< 2 



whenever < \x\ < r\. Let R 2 = T\ A L 1 1 . Then for x € M d \ {0} we have 
G(x)-G(Lix) > 



\x\ d - a e(\x\- 2 ) \L lX \ d - a £(\L lX \- 2 ) 
c 2 £(\x\- 2 ) 



1 



> 



\x\ d - a £(\x\- 2 ) V L d - a e(\L!x\- 2 ) 
1 1 



2c \x\ d - a £{\x\~ 2 ) 



n 



Proposition 3.2 for every bounded open set D, there exists a constant C\q > 1 such that 

1 



G D {x,y) > C^ 1 



G D (x,y)<C w 



1 



for x,y e D . 



(3.1) 



\x - y\ d - a £(\x ~ y\~ 2 ) ' 

for x, y £ D with Li\x — y\ < Sd(x) A 5o(y), (3-2) 



\x-y\ d -«l{\x-y\- 2 y 
where L\ is the constant in Lemma \3.1l 

Proof. Since Grj(x,y) < G(x,y), D is bounded and £ locally bounded, (|3.ip is an immediate consequence 
of Theorem 12.91 Now we show ()3.2|) . Without loss of generality, we assume 8d{v) < o~d{x), and let 
M := diam(Z3). We consider three cases separately: 

(a) 8 D (y) < R 2 : Since \x-y\ < 8 D (y)/L l7 \X TBiy5D{y)) - y\ > S D {y) > L x \x- y\. Thus by the monotonicity 
of G and Lemma EUl 



G D {x,y) > G B{y ,s D ( y) ){x,y) = G(x,y)-E x \G{X TB{ySD{y)V y) 

1 



> G{x - y) - G{L x {x - y)) > c x 



\x - y\ d ~ a £(\x - y\~ 2 )' 



V|a;-y| < 



5d(v) 



(b) Sn(y) > R 2 and \x — y\ < R<ijL\: In this case, \X Tb{ r —y\ > R 2 > L\\x — y\ and, by the monotonicity 
of G and Lemma 13. 1[ we get 



G D (x,y) > G B{v<R2) (x,y) = G(x,y)-E x \G{X TB(yIi2V y) 

1 



> G(x-y)-G(L 1 (x-y)) > a 



R 2 



\x-y\ d -«£{\x-y\- 2 Y ^ yl ~ L x 



(c) 8n{y) > R 2 and \x — y\ > R2/L1: Choose a point w £ dB(y,R2/Li). Then from the argument in (b), 
we get 



G D {w,y) > G(w,y)-E w G{X TB{yR2) ,y) 



> ci- 



1 



(R 2 /L 1 ) d -£((R 2 /L 1 )- 2 )- 

Since D is bounded and Grj(- ,y) is harmonic with respect to X in B(x,R 2 /(2L 1 )) U B(w,R 2 /(2Li)), by 
Theorem 12. 141 we have 



Go(x,y) > c 2 G D {w,y) > c 3 



1 



{R 2 /L 1 ) d -^£{{R 2 /L 1 )- 2 ) 
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* k(rX)->) L<i<M i{s ~ 2) ) ix-»i--4ix-»i- a )' v|s - y|> f- 



□ 



Lemma 3.3 for every L > and bounded open set D, there exists C\\ > suc/i f/iat /or every \x — y\ < 
L(6 D (x)A6 D (v)), 

Proof. Without loss of generality, we assume 5/j(x) < 6d{v)- Moreover, by Proposition 13.21 we can assume 
that L > 1/Li and we only need to show (|3.3j) for j-(5u(x) < \x — y\ < LSd(x). 
Choose a point w £ dB{x,8D(x)/L\), Then by Proposition [321 we get 

Gd(x,w) > ci- 



{S D (x)/L 1 y-"i((S D {x)/L 1 )- 3 )' 

Since |y — u;| < \x — y\ + \x — w\ < (L+ 1)8d(x) and Gd(x, • ) = Gd( ■ >#) is harmonic with respect to X in 
B(y,S D (x)/Lx) UB(w,S D (x)/Li), by Theorem l2~T4l we have 

i 
Gd{x,v) > c 2 G D (x,w) > c 3 



(6 D (x)/L 1 y-^£((5 D (x)/L 1 )~ 2 ) 



- C4 U(M*)/£i)- 2 ) J |a; - W^iQx - y\~ 2 ) ' (3-4) 

By the uniform convergence theorem ([2j Theorem 1.2.1]), we can choose a small ri = ri(£, L) > such that 

bf «££>>I. W< ri . (3.5) 

Ae[i,2i] l(r~ 2 ) 2 

If ±6 D (x) <\x-y\< LS D (x) < n, by @3M£5), 

rr ^> f inf i(W 2 A x ^ c 4 i 

Dl ' y; - 4 Vr<n,A€[i,2i] e(r- 2 ) J \x - y\ d - a £(\x - y\~ 2 ) - 2 \x - y\ d ~ a e(\x - y\- 2 )' 

On the other hand, if -j-8d{x) < \x — y\ < LSd(x) and Srj(x) > ri/L\, then |aj — J/| > -r^ - - Thus from 
(|3.4p . we see that 

G fl (i,i/)>c 4 inf £(r- 2 ) j ( inf ^(r^" 1 ' 



^el-^.A/] y \r&[^,M] J \x-y\ d - a l{\x-y\~ 2 ) 

where M — diam(D). D 

For the remainder of this section, we assume that D is a bounded K-fat open set with characteristics 
(Ri,k). Without loss of generality we may assume that i?i < 1/4. We recall that for each Q £ dD and 
r £ (0, Ri), A r (Q) denotes a point in D n B(Q, r) satisfying B(A r (Q), nr) C D n B(Q, r). Recall also that 
Gd{-, z) is regular harmonic in D \ B{y, e) for every e > and vanishes outside D. From Theorem l2.151 we 
get the following boundary Harnack principle for the Green function of X in D. 

Theorem 3.4 There exists a constant Cn > 1 such that for any Q £ dD, r £ (0, R%] and z,w £ D\ 
B(Q, 2r), we have 

! G D {z,A r (Q)) G D {z,x) G D {z,A r {Q)) cnn n(nL\ 

° 12 G D (w,A r (Q))- G D (w,x) -^ 12 G D {w,A r {Q)Y ^^ n ^> 2 J ■ 
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Using the uniform convergence theorem ([U Theorem 1.2.1]), we further choose R3 < R\ such that 

'<_ mi „ «<^i< imaXi ^)< 2 , itr <« 3 . (3.6, 

Let us fix zq £ D with kR^ < Sd(zo) < R3 and let E\ :— nR^/24,. For x,y £ D, we define r(x,y) := 
6d(x) V 5 D (y) V\x — y\ and 

{{AeD: S d (A) > %r(x,y), \x - A\ V \y - A\ <5r{x,y)} iir{x,y)<e 1 
\{z Q } if r(x,y) > e\. 

Note that for every (x, y) £ D x D with r(x, y) < e\ 

Ud(A) < 5 D {x) V 6 D (y) V\x-y\ < 2 K - 1 5 D (A), A £ B(x, y). (3.8) 
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So by flSU), ifr(»,y) < £1 

1 £(fe(A)) - 2 ) 

2 " : ^(rfoy)) 

Let 



^ 777I7I n^3^: < 2 > A£B(x,y). (3.9) 



C* 13 :=C 1 o2 d - Q ( 5 D (z )- d +« sup ^r" 2 )" 1 . (3.10) 

\M*o)/2<r<Af y 

It follows from Proposition 13.21 that Gd(-,zq) is bounded above by C13 on D \ B(zq, Sd(zo)/2). Now we 
define 

5 (a;):=G D (x,z )ACi3. (3.11) 

Note that if #0(2) < 6ei, then \z — z \ > 8b{zq) — &£\ > Sd(zq)/2 since 6ei < 6d(zq)/4, and therefore 
#(z) = G D (z,z Q ). 

The two lemmas below follow immediately from Theorem 12.141 

Lemma 3.5 There exists C14 = Cu(k) > 1 such that for all x,y £ D, C^g{A\) < 5(^2) < Cng{Ax) for 
all Ai,A 2 £ B(x,y). 

Lemma 3.6 There exists C15 > 1 such that for every x £ {y £ D : 5jj{y) > k 3 £ 1 /64}, C^ 1 < g(x) < C15. 



We are now ready to prove Theorem II .21 

Proof of Theorem 11.21 Since the proof is an adaptation of the proofs of 4, Proposition 6] and [TBI 

Theorem 2.4], we only give the proof when <5_d(cc) < #o(y) < f \x — y\. 

In this case, we have r(x, y) = \x — y\. Let r := |(|a; — j/| Aei). Choose Q x , Q y £ dD with \Q X ~ x\ = 8d{x) 

and \Q y — y\ = fo(y). Pick points xi = A Kr / 2 (Q x ) and 3/1 = A rer / 2 (Qy) so that x,xi e B(Q x ,nr/2) and 

y, j/i G B(Q y , nr/2). Then one can easily check that \z — Q x \ > r and \y — Q x \ > r. So by Theorem 13.41 we 

have 

_ x G D (xi,y) G D {x,y) G D {x l7 y) 

g(xi) g{x) g[xi) 

for some c\ > 1. On the other hand, since \zq — Q y \ > r and \x\ — Q y \ > r, applying Theorem 13.41 again. 

-1 Gpfayyi) < G D (x 1 ,y) Gp(ii,»/i) 

1 g(yi) g(y) ~ 1 3(2/1) 
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Putting the two inequalities above together we get 

c _ 2 G D (xi,yi) G D (x,y) 2 G D (x 1 ,y 1 ) 
1 g{xi)g{yi) ~ g(x)g(y) ~ 1 g(x 1 )g(y 1 )' 

Moreover, ||a; — y\ < \x x — yi\ < 2\x — y\ and \x x — yi\ < -^t(5d(x X ) A 5n{yi))- Thus by Lemma |3"UI we 
have 

2- d +«c 2 - 1 cr 2 1 < G D {x,y) < 3 d -°c 2 cf 1 



g{xi)g{yx) \x - y\ d - a £(\x x - y x \~ 2 ) " g(x)g(y) ' g(x x )g(y x ) \x - y\ d a t(\x x -yi\ 2 ) 

for some c 2 > 1. 

If r = £i/2, then r(x,y) = \x — y\ > £ X . Thus g(A) = g(zo) — C13 and 5d{xi) A Soiyi) > Kr/2 = ne x /A 
So by Lemma l3.6[ 

C1-3V < -r^n-, < ^iV 2 (3.13) 



for some C3 > 1. 

If r < £i/2, then r(x,y) = |ar — y| < £1 and r = |r(o;, y). Hence £0(21), fo (2/1) > Kr/2 = nr(x,y)/4. 
Moreover, \xi — A\, \yi — A\ > 6r(x, y). So by applying Theorem 12.131 to g twice, 

c -i < 44 < c 4 and c- 4 x < 44 < c 4 (3.14) 

for some constant C4 = 04(D) > 0. Combining (|3.12|l - (|3.14|l . we get 

-1 g(x)gjy) r , , g(x)g{y) , R , s 

Cs g(AY\x - y\ d ~H{\x x - yih 2 ) " D[X ' V) ~ ° 5 g(A)*\x - y\ d -<*e(\x x - y x \-*) ' A & ° {X ' y) ' 

If §|as - y| < In - 2/i| < 2|.x -y\< R 3 , by (JSU), § < f^Z^-l) < 2- On the other hand, if \\x - y\ < 
\xi -yi\<2\x- y\ and 2\x -y\> R3, we have that -^ < \x - y\ < M and &■ < \xi - y x \ < M. Thus from 
the local boundedness of t on (0, 00), we see that 

-1 < l(|s-y|- 2 ) < 

Ce - e(\x x - y2 \- 2 ) - C6 

for some eg > 1. □ 



4 Explicit Green function estimates on bounded C ll -open sets 

In this section we refine the estimates from Theorem 11.21 in the case of bounded C 1,1 open sets. 

Recall that X — (Xt : t > 0) is the d-dimensional subordinate Brownian motion defined by Xt — Ws t 
where W = (W 1 , . . . , W d ) is a d-dimensional Brownian motion and S = (St : t > 0) an independent 
subordinator with the Laplace exponent which is a complete Bernstein function satisfying assumption 
(H) for a e (0, 2 A d) and the additional assumption (|2.14|) when d < 2. Let Z = (Z t : t > 0) be the 
one-dimensional subordinate Brownian motion defined as Z t :— Wg . Recall that the potential measure of 
the ladder height process of Z is denoted by V and its density by v. We also use V to denote the renewal 
function of the ladder height process of Z. In Corollary |2.8l we have established that both V and v are C°° 
function. Recall the notation R^ := {x = (x x , . . . ,Xd- X ,Xd) :— (x,Xd) G f 1 : ^ > 0} for the half-space. 
The next result, which follows from |28[ Theorem 2], is very important in this paper. 

Define w(x) := V((x d )+). 
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Theorem 4.1 The function w is harmonic in K.1 with respect to X and, for any r > 0, regular harmonic 



pd-l 



x (0,r) fm (0,r) w/ien d — 1) with respect to X . 



Proof. Since Z t := W$ has a transition density, it satisfies the condition ACC in [28], namely the resolvent 
kernels are absolutely continuous. The assumption in 28 that is regular for (0, oo) is also satished since Z 
is symmetric and has infinite Levy measure. Indeed, if were irregular for (0, oo), it would be, by symmetry, 
irregular for (— oo, 0) as well. But then Z would be a compound Poisson process which contradicts the 
fact that it has infinite Levy measure. Further, again by symmetry of Z, the notions of coharmonic and 
harmonic functions coincide. Let Z^ 0,00 ) (respectively X R +) denote the process Z killed upon exiting (0, oo) 
(respectively X killed upon exiting M.+ ). By [28j Theorem 2], the renewal function V of the ladder height 
process of Z is invariant for Z^- '°°> . Thus w is invariant for X + . Being invariant for X + , w is also harmonic 
for X R +, and consequently, harmonic in R+ with respect to X. We show now that w is regular harmonic for 
X in R d_1 x (0, r) for any r > 0. First note that since V is continuous at zero and 1^(0) = 0, it follows that 



lim w(x) — lim w(x,Xd) — lim V{xd) =0. 

IEj-S-0 zj-s-O Xrf-i-0 



(4.1) 



Thus, by harmonicity of w and (|4.1 



w(x) = w{x,x d ) = limE x w \X TRdl ^^ \ = E x w ( 



^"'xlO.r) 



x d > . 



□ 



Proposition 4.2 For all positive constants rg and L, we have 



sup / w{y)j{\x-y\)dy < oo . 

xGR d : 0<x d <L J B(x,r ) c nR^ 

Proof. By Theorem |4~T1 and (|2.19j) . for every jieRl 



w(x) > E a 



w{X T 



'('■i)/2)n«: 



)'X T , 



(*,r /2)nE: 



d eB(x,r ) c r)R d + 



B(x,r ) c r\R± JB(x,ro/2)nRi 



G B(x,r /2)nR d (x, z)j(\z - y\)w(y) dz dy . 



(4.2) 



Since \z — y\ < \x— z\ + \x — y\ < r + \x — y\ < 2\x — y\ for (z,y) e B(x,r /2) x B(x, r ) c , using f|2.17|) and 
(|2.18j) . we have j(\z — y\) > c\j{\x — y\). Thus, combining this with (|4.2|) . we obtain that 

f t \ •/! IN I -1 W(X) 

sup / w(y)j(\x~y\)dy < Cl sup — =-. 

x£R d :0<x d <LJ S(s,r ) c nR^ x=0,0<x d <L ^x [ T B(x,r /2)nR d _ \ 

We claim that the supremum on the right-hand side above is finite. Clearly, if L > Xd > ro/(64) and J = 0, 

w(x) < y(£) 

E cc[Ts( 2; ,r /2)nK^] ~~ E [TB(0,t- o /(64))] 

Suppose Xd < r /(64) and x = 0. Let U :— B((0, 16r ), r ). By the Levy system, we have 



11 x [Xr d e U 

1 B(0,r /2)nlt» 



(7 JB(0,r /2 



G i3(0,ro/2)nRi( ;c . z )j(l 2: - 3/|)d«d» < C2E x [T B(0tro/2)nM d 
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Thus, by the above and the boundary Harnack principle (Theorem 12. 15|) . 

w(x) w(x) w(xi) Tr , ,, „ ss 

< C 2^T^ — TTn < c 3^^ — "Fn < c*V(ro/(16)) 



^x\T R(7 . rn/2 ) nR d} ~ P X (X T . €U) - F X JX T „ eU) 

where x\ = (0, ro/(16)). We have thus proved the claim. □ 

We now define the operator (A, T)(A)) by the following formula: 
Af(x) := P.V. I (f(y)-f(x))j(\y-x\)dy:=lim[ (f(y) - f(x)) j(\y - x\) dy 

JR d £ ^' J{y£R d :\x-y\>e} 



D(A) := < / :R d -»R:lim / (f(y) - f(x))j(\y - x\) dy exists and is finite > . (4.3) 

[ £ i° J{ v m d :\x-y\>e} J 

It is well known that Cq C £>(_4), where Cq is the collection of C 2 functions in R d vanishing at infinity, and 
that by the rotational symmetry of X, A restricted to Cq coincides with the infinitesimal generator of the 
process X (e.g. [26l Theorem 31.5]). 

Theorem 4.3 Aw(x) is well defined and Aw(x) = for all x £ R+. 



Proof. We first note that it follows from Proposition 14. 21 and the fact that j is a Levy density that for any 
£>0 



/ (w(y)-w(x))j(\y-x\)dy 



sup 

xef: 0<x d <L 

< sup / w(y)j(\y-x\)dy + V(L) I j(\y\)dy <oo. (4.4) 

i6l d : 0<x d <L J{y<ER d :\y-x\>l} J {y£R d :\y\ >1} 

Hence, for every e € (0, 1/2) 

„4 £ w(x) := / (w(y) - w(x))j(\y - x\)dy 

J {y£R d :\y-x\>e} 

is well defined. Note that since w(x) — V((xd) + ) and V is smooth in (0, oo) by Corollary 12.81 it holds that 
w is smooth in Mi_. Hence, 



A s w(x) = / (w(y) - w{x) - l{\ v - x \<\}{y - x) • Vw(x)) j(\y - x\)dy. 

J {y£R d :\y~x\>e} 



Moreover, by the smoothness of w, 



X M- / (w(y) - w(x) -(y-x)- Vw(x)) j(\y - x\) dy 

J {yGR d :\y-x\<e} 

converges to locally uniformly in R+ as e — > . Combining this with (|4.4j) . we see that Aw is well defined 
in R+ and A E w(x) converges to 

Aw{x) = / (w(y) - w{x) - l{| J/ _x|<i}(2/ - x) ■ Vw(x)) j(\y - x\)dy 

JR d 

locally uniformly in R^ as e — > 0. 
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Moreover, for every x £ R+, z £ B(x, (eAXd)/2), and y £ B(z,e) c it holds that h\y— z\ < \y— x\ < ||y— z| 
So, using ([2~T7j) . 

l{|y- z |> £ } |0(y) - u>(z) - l { | y _ z |<i } (y - z) ■ Vw(z))\ j(\y - z\) 
<c\ sup V"(a) )\y-x\ 2 l {e / 2< \y^ x \ < 2}j{\y-x\/2) 

\e/2<s<x d +2 ) 

+ (w(y) + V(x d + l))l {ly - xl>1/2} j(\y - x\/2). 



It follows from Proposition 14.21 and the fact that j is a Levy density, by using the dominated convergence 
theorem, that x — > A e w(x) is continuous for each e. Therefore, by this and the local uniform convergence of 
A e w, the function Aw(x) is continuous in W.i_. 

Suppose that U\ and Ui are relatively compact open subsets of R^ such that U\ C U2 C U2 C R+. Let 
r := dist(C/i, [/J) > 0. Then, by Proposition [ 



^(y)i(k _ y\)dydx < \Ui\ sup / w(y)j(|a: - j/|)dj/ 
UxJus xEUiJm 



< |E7i| sup / w(y)j(|x - yDdy < 00. (4.5) 

xGE/i JB(x,r )<= 

By harmonicity of u>, clearly w(X TU ) <G L 1 (P X ) and 

sup E^ [l u: ={X TUi )w(X TUi )\ < sup E x [w(X TUi )] = sup iu(x) < 00. 

The last two displays show that the conditions [7] (2.4), (2.6)] are true. Thus, by [3 Lemma 2.3, Theorem 
2.11(h)], we have that for any / £ C 2 (R^), 

0=/ f (w(y)~w(x))(f(y)-f(x))j(\y-x\)dxdy. (4.6) 

JR d JR d 

For / £ C C 2 (R^) with supp(/) C U[ C U 2 C % C R f { , 

|u>(y) - u>(x)||/(y) - /(a;)|j(|3/ - a;|)efod£v 



|w(y) - w{x)\\f{y) - f(x)\j(\y-x\)dxdy + 2 / / \w(y) - w(x)\\f(x)\j(\y - x\)dxdy 
U 2 JU2 JUi Ju% 

<ci / |y-x| 2 j(|y-a;|)da;dy H- 2 1 1 / 1 1 oo | t/ x | ( sup w(x) ) / j{\y - x\)dy 

JU2XU2 \xeUi J JUS 



+ 2||/||oo / / ^(2/)i(|x-y|)dy& 

is finite by (|4.5|) and the fact that j is a Levy density. Thus by (|4.6j) . Fubini's theorem and the dominated 
convergence theorem, we have for any / £ C 2 (Ml ) , 



= lim / (w(y) -w(x))(f(y) - f(x))j(\y - x\)dxdy 

£ ^° J{(x,y)eR d xR d , \y-x\>s} 

= -21im/ f{x)[ / (w(y) - w(x))j(\y - x\)dy ) dx = -2 / /(a;)^4u;(aj) cte, 

where we have used the fact A e w — > Aw converges uniformly on the support of /. Hence, by the continuity 
of Aw, we have Aw(x) — in R^. □ 
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For x £ R d , let Sgoix) denote the Euclidean distance between x and dD. It is well known that any C 1,1 
open set D with characteristics (R, A) satisfies both the uniform interior ball condition and the uniform 
exterior ball condition with the radius r%: there exists r% < R such that for every x £ D with Sjj{x) < r\ 
and y £ M. d \ D with 5n(y) < ri, there are z x , z y £ dD so that |x — z x \ = 8qd{%), \y ~ z y \ = 6d(.v) and that 
S(x ,ri) C D and B{y ,r t ) cR d \D for x = z x + r x {x- Z x )/\x- z x \ and y Q = z y + r x {y- z y )/\y- z y \. 

In the remainder of this section, we assume D is a C 1,1 open set with characteristics (i?, A) and D satisfies 
the uniform interior ball condition and the uniform exterior ball condition with the radius R (by choosing 
R smaller if necessary) . 

Lemma 4.4 Fix Q £ dD and let 

h(y) := V (S D (y)) lDnB(Q,R)(y)- 

There exist C\q — Cie(a, A, R,£) > and R4 < R/A independent of the point Q £ dD such that Ah is well 
defined in D D B{Q, R4) and 

\Ah{x)\<C w for allx£ DnB(Q,R±). (4.7) 



Proof. We first note that when d = 1, the lemma follows from Proposition 14.21 and Theorem 14.31 In 
fact, suppose that d = 1 and x £ D n B(Q,R/2). Without loss of generality we may assume that Q 
is the origin and D PI B(Q,R) = (0,R) (due to uniform exterior ball condition). Since h(y) — w(y) for 
y £ D n B{Q, R) = (0, R), we have 



A{h-w){x) = lim/ (^ — tw)(y)i(|y — x|)dj/ 

eJ-0 7{v6H I :|a:-v|>e} 



lim / (h-w)(y)j{\y-x\)dy 

£i0 J{y&(0,R)":\x-y\>e} 



lim / w(y)j(\y-x\)dy 

£ i° J{yG(0,fl) c :|x-y|>e} 



= -lip / w(y)j(\y-x\)dy 

e-l-0 J {y>R:\x-y\>e} 

Since < a; < i?/2 and y > R, we have |x — y| > -R/2, and thus 

\A(h-w)(x)\ < / w(y)j(|y-2:|)ck/. 

^{y> J R,|x- J/ |> J R/2} 

Therefore, by using Theorem l4.3l (which gives ^4/i(x) = A{h — w)(x)), Proposition ^. 21 and the above display, 
we conclude that 



\Ah(x)\ = \A(h - w)(x)\ < sup / w{y)j{\y- z\)dy < 00. 

0<z<.R/2 J{y£R^:|y-z|>-R/2} 

Throughout the remainder of the proof, d > 2 and q is a fixed positive constant such that 

0<q< aA (2 ~ a) and a + 2<z-1^0. 

Since £ is slowly varying at 00, by Potter's Theorem ([H Theorem 1.5.6 (1)]), we can find a small 
Ri < 1 A {R/A) such that for every r < 2R\ 

^ 2 ) < 2 £ (( 2J? D~ 2 ) /OR2X1/2 -1/2 < -1/2 f48) 



(£((2i?- 1 )- 2 r~ 4 )) 1 /2 - (^((2i?- 1 )- 2 (2i?l)- 4 )) 1 / 
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1{M) < 2£((2R 2 4 )- 1 )(2Rlyr- q <Cir- q , (4.9) 

t{r~ 2 )- 1/2 < 2i{{2Rl)- 2 ) 1 ' 2 {2Rlfr- (i < Cl r- q . (4.10) 

In the remainder of this proof, we fix x 6 D n B{Q,R/C) and xq G dD satisfying 8d{x) = \x — Xq\. We 
also fix the C 1 ' 1 function ip and the coordinate system CS = CS Xo in the definition of C 1,1 open set so that 
.x = (0, x d ) with < x d < R 4 and B(x , R) D D = {y = (y, y d ) e B(0, R) in CS :y d > i/j(y)}. Let 



Mv)-=R-\/R 2 -\y\ 2 and V 2 (y):=-i?+^ 2 -|yl 2 - 

Due to the uniform interior ball condition and the uniform exterior ball condition with the radius R, we have 

MV) < ^{y) < My) for every y£ DnB(x, i? 4 )- (4.11) 

Define H + := {y = (y, y d ) £ CS : y d > 0} and let 

A := {y = (y, y d ) e (D U H+) n B(x, R A ) :i> 2 (y)<y d < ^i(y)}, 

^ := {y = {y,Vd) e b{x,r a ) ■. yd > My)}- 

Note that, since |y — Q| < |y — x| + \x — Q\ < R/2 for y 6 i?(a;, i?4), we have 

B(x,Ri)nD cB{Q,R/2)nD. (4.12) 

Let 

Note that /i x (x) = h(x). Moreover, since S , (y) — (y d ) + in CS, by Theorem 14.31 it follows that Ah x is well 
defined in H + and 

Ah x (y) = 0, VyeH+. (4.13) 

We show now that A(h — h x )(x) is well defined. For each e > we have that 

(h(y) - h x (y))j(\y - x\) dy 

{y£DUH+: | y — a;|>e} 

< / (Hy) + h x {y))j{\y - x\)dy + f {h{y) + h x (y))j{\y - x\) dy + f \h(y) - h x (y)\j{\y - x\)dy 

JB(x,Ri)" J A JE 

=:h+h+h- 

We claim that 

h+h+h<C w (4.14) 

for some constant C\q = Cxe(a, A, R,£). This shows in particular that the limit 



lim / (h(y) - h x (y))j(\y - x\) dy 

e ^' J{y£DuH+:\y-x\>e} 

exists and hence A(h — h x )(x) is well defined, and \A(h — h x )(x)\ < C\q. By linearity and ()4.13|) . we get 
that Ah(x) is well defined and |.A/i(x)| < Cx§. Therefore, it remains to prove (|4.14|l . 
By the fact that h(y) = for y e B(Q,R) C , 

h< sup / V(y d )j(\z-y\)dy + c 3 j(\y\)dy =: K x + K 2 . 

2£R d : Q<z d <RJ B(z,Ri) c nH+ J B(0,_R 4 ) c 

Ki is clearly finite since J is the Levy density of X and K\ is finite by Proposition ^. 21 
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For y £ A, since V is increasing and (R — \JR 2 — |y| 2 ) < R 1 \y\ 2 , we see that 

h x (y) + h(y) < 2V(Mv) ~ MV)) < 2U(2iT 1 |y| 2 ). (4.15) 

Using (|4TT5|) . (|2~T2"T) and Theorem I2TH we have 

h < I * / lA(y)(fcx(y) + %M(r 2 + |j/rf - a: d | 2 ) 1/2 ) md -i(dy)dr 
JO */|y|=r 

/•-R4 /• 

<2 / / l A (2/)F(2i?-V 2 )j(r) m d -i(dy)dr 

JO J\y\=r 

<^J o r ~ d m J R 2 r Ly/2 m *-i({y eA'-m= r})dr (4.16) 

where rrid-i is the surface measure, that is, the (d — l)-dimcnsional Lebesgue measure. Furthermore, since 
|^2(y) - "01 (y ) I < 2R- l \y\ 2 = 2R- l r 2 on \y\ = r, we have for r < i? 4 , m d _!({y : |y| = r,ip 2 (y) < y d < 
V'i(y)}) < c 5 rd f° r some constant C5. Using the above inequality and (|4.8JI . from (|4.16|) we get 

/■H4 ^(r -2 ) f^ 4 

/2 - C6 io W(2fl-i)-^))i/ a dr ^i, r- 1/2 *«»- 

For ^3, we consider two cases separately: If < y c i — S + (y) < 5d(v), since i> is decreasing, 

%) - K(y) < V(y d + R-^y] 2 ) - V(y d ) = / v(z)dz < R- l \y\ 2 v{y d ). (4.17) 

Jyd 

If yd = 6 , (y) > 8d{v) and y E E, using the fact that &d{v) is greater than or equal to the distance between 
y and the graph of -01 and 



y d -R + V\W + (R-y d ) 2 = ^== =-^ lu -<X7^L<^ 

^|y| 2 + (it!-y d ) 2 + (i?-y d ) 2{R-y d ) R 

we have 

yd 

1 1 t.\ 2 , 



h x (y)-h(y)< v (z)ciz < iT ^l 2 v(R - y/\y\ 2 + (R- y d ) 2 ). (4.18) 

J R-yJ\W+(R-y d ) 2 



By dum-dsg), 



h<R~ X \y\ 2 v{y d )j{\x-y\)dy 

J En{y.y d <S D (y)} 



+ R- 1 / \y\MR - V\W + (R-yd) 2 )j(\x - y\)dy =■. R-\L X + l 2 ). 

J En{y.y d >S D (y)} 

Since E C {z = (z,z d ) e R d : \z\ < R 4 and < z d < 2i? 4 }, using polar coordinates for y and by 
Theorem [2~TT1 (|2~T2l) . (jO|) and CTO)) . we have that 

,-2i?4 / /--R.4 \ 

ii < <*j v(y d )(j r 2 j{{r 2 + \y d -x d \ 2 ) 1 ' 2 )r d - 2 dr\dy d 

f 2 Ri 1 ( [*< r^((r 2 + |y d -.T rf | 2 )- 1 ) J \ 



(yd) 1 -"/ 2 ^^ 2 )) 1 / 2 Wo (^ 2 + l^-^l 2 ) (d+a)/2 
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2fl 4 i / r R i 



~ CW Jo (yd) 1 -^ 2 ^ IX (r2 + \y d -x d \2)(^+^ dr ) dyd 



2/l4 -1 / /"^14 1 



- "°/o (v,) 1 -/ 9 * U FTk^TF^'V dyd 

2fl 4 1/1 1 



- Cn X (yd) 1 -"/ 2 ^ V |W " ardl^ 2 *- 1 + (H4 + |Vd - x d |)«+ 2 9-i ) dyd ~ ° 12 

for some constant c&, . . . , c\i > 0. The last inequality is due to the fact that q < (2 — a)/20, which implies 
(1 — a/2) + a + 3q — 1 < (6 + 7a)/20 < 1, so by the dominated convergence theorem, 

/•2R4 -L 

^X (^W^-s^-i^ (4A9) 

is a strictly positive continuous function in x d G [0, R4] and hence it is bounded. 

On the other hand, we have, using polar coordinates for y, and by Theorem 12. 11[ (|2.12[) and (|4.9|) - (|4.10j) . 



x d +R A / rR4/\y/2Ry d -y 2 d 



L 2 <c 13 J I / v(R - Vr 2 + (R- Vd) 2 )r d j((r 2 + \y d - x d \ 2 f' 2 )dr dy d 



( z d +R, J rIUKy/2Ry d -v> {R _ ^ r 2 + {R _ yd) 2 )a /2- H{{r 2 + | _ ^2)-!) \ 

< C14 / / ; ■ - — —r dr dyd 

Jo \Jo (i((R - V r2 + (R- yd) 2 )- 2 ))^ 2 (r 2 + \y d - x d \ 2 )( d +«)/ 2 J 

j-x d +Ri / i.R i /\ y /2Ry d -y 2 d ^d \ 

Jo V° ( R - V r2 + {R~ yd) 2 ) 1 -*' 2 ^ (r 2 + \y d - Xd \^)(d+a+2 q )/2 J 

i-x d +Ri I /•Ri/\y / 2Ry d -y 2 l ^ \ 

" Cw Jo IX (R - V 2 + (R ^yd?Y' a,2+q (r + \vd - xd\) a+2 o d J dVd ' 



Since, for < r < R 4 A y/2Ry d - y 2 ,, 



i _ R + v 7 ^ 2 + (R- Vd) 2 K ci 7 



i? - ^ + (i? - 2/d) 2 W2Ry d -yj + r)(^/2Ry d -yj-r) " y/m(y/2Ry d - y 2 d - r) ' 
we have 



l-Xd+Ri Cig rRi^2Ry d -y 2 d ^ 

L2 - X (2/ d ) (1 - Q/2+(z)/2 X ( V 2 ^ - ^ - r) W2+ ? (r + |y d - a: d |)<*+ 2 9 ^ ' 



Using the fact that q < ^, we see that with a := \j2Ry d — y 2 d and b := |y d 



»d|, 



1, 



i?4Aa 



rfr 



(a - r ) 1 -«/ 2 +9 (r + b) a + 2 i 

(i?4Aa)/2 dr ,i? 4 Aa ^ 



(a - r) 1 -"/^ (r + b) a + 2 i J {RiAa)/2 (a - r) l - a / 2 +<i (r + b) a + 2 <i 

2l-a/2+q f(Rif\a)/2 ^ ^ ^R A r\a ^ 



a 



l-a/2+q 







(r + b)«+ 2 <i (b + (i? 4 A a)/2)-+ 2 i 7 (%Aa)/2 (a - r)i-«/ 2 +? 



<- £19 , C19 a/2-g <- £20 C 2 Q 

- a l- Q /2+ ?& (a+2q-l)+ + (^ A a )a+2 9 - ( yrf ) (l-a/2+,)/2 \ Xd _ y<J | (a+2,-l)+ + ( yd )(a+6 ? )/4 
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Thus we obtain 

r2Ri dy d f 2Ri dy d 

{yd)( i- a /2 +q) \ yd _ Xd \( a+ 2 q -i) + + c 2i J q (^)(I+4 g )/ 



L 2 < 021 / ~, NM _„ /24 .„i|., „ U n +2n-U+ +C21 / /■■ UU^/2 - ( 4 - 20 ) 



Since q < 1/10, the second integral in (J4.20I) is bounded. And by the same argument as the one for (|4.19p . 
the first integral in (|4.20p is also bounded. We have proved the claim (J4.14I) . □ 

When d > 2, define Pq{x) := x d — iPq(x), where (x, x d ) are the coordinates of x in CSq. Note that for 
every Q £ dD and x £ B(Q, R) n D we have 

(l + A 2 )- 1/2 PQ(x) < S D (x) < PQ {x). (4.21) 

For a, b > 0, we define Dq(o, b) := {y £ D : a > pcz{y) > 0, \y\ < b} when d > 2. When d = 1, we simply 
take D Q (a, b) = D Q (a) := B(Q, a) n D. 



Lemma 4.5 T/iere are constants R 5 = R 5 (R, A,a,£) £ (0, i? 4 /(4 v /l + (1 + A) 2 )) and d = C;(5,A, a) > 0, 
i = 17, 18, such that for every r < R§, Q £ dD and x £ Dq{t, r), 



J x (x TDQ(rr} £D)> C 17 V(5 D (x)) (4.22) 



and 

E, [r DQ (r,r)] < C ls V(S D (x)). (4.23) 

Proof. Without loss of generality, we assume Q = 0. For <i > 2, let ^ = ^ : R d_1 — > R be the C 1 * 1 function 
and CSq be the coordinate system in the definition of C 1 ' 1 open set so that 5(0,5) (ID = Uy, y d ) £ 
5(0,5) in CS'o : y d > ^(y)}. Let p(y) := y d - ^(£) and D(a,6) := D (o,6) for d > 2. When d = 1, D(a,6) 
is simply B(Q, a) n D. The remainder of the proof is written for d > 2. The interpretation in the case d = 1 
is obvious. 
Note that 

\V\ 2 = \V\ 2 + \Vd\ 2 <r 2 + (\y d - i>(y)\ + \^(y)\) 2 <(! + (! + A) 2 )r 2 for every y £ D(r, r) . (4.24) 

Hence, by letting f := R i / y /l + (1 + A) 2 , D(r, s) C D(r, r) C 5(0, 5 4 ) n 5 C 5(0, 5) n 5 for every r,s <r. 
Define 

%) : = V(6 D (y))iB(o,R)r\D(y)- 

Let <? be a non-negative smooth radial function with compact support such that g(x) = for \x\ > 1 and 
J Rd g(x)dx = 1. For fc > 1, define ^(z) = 2 kd g(2 k x) and 

/i (fc) (z) := (g k * h)(«) := / g k {y)h{z - y)dy , 



and let 5 fc := {a; e D n 5(0, 5 4 ) : <Wj(o„r 4 )(z) > 2~ fc }. Since h (fc ) is C°°, „4M fc ) is well defined everywhere. 
We claim that 

- C 16 < Ah,™ < C w on B k , (4.25) 

where C\$ is the constant from Lemma [4.41 Indeed, for x £ B k and z £ 5(0, 2~ fc ) it holds that x — z £ 
D n 5(0, 5 4 ). Hence, by Lemma \A. 41 the following limit exists: 



e ^° -'\x-y\>e 



lim / (h(y-z)- h(x - z)) j(\x-y\)dy 
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lim 



(h(y') - h(x - z)) j(\(x - z) - y'\) dy' = Ah(x - z) 



>\{x-z)-y'\>e 

Moreover, by the same Lemma B~4l it holds that — C 16 < Ah(x — z) < C w . Next, 



(hW(y)-hW(x))j(\x-y\)dy 



x—y\>e 



\x—y\>e 



9k(z)(h(y - z) - h(x - z)) dz j(\x - y\) dy 



\z\<2- 



9k{z) 



(h(y - z) - h(x - z)) j(\x - 3/|) dy dz. 



x-y\>e 



By letting e — > and using the dominated convergence theorem, it follows that 

Ah ( - k \x) = / g k (z)Ah(x - z) dz < C 16 / g k (z) dz = C 16 

J\z\<2- k J\z\<2- k 



The left-hand side inequality in (|4.25j) is obtained in the same way. 

Using the fact that A restricted to C^° coincides with the infinitesimal generator of the process X, we 
see that the following Dynkin formula is true; for / <E C^°(M ) and any bounded open subset U of R d , 



E, / Af(X t )dt = E x [f(X ru )} - f{x). 



(4.26) 



Let U C D n B(0, R 4 ). By using (Ei~2"oT) for U D B k and h^ k \ the estimates P~2"5j) . the fact that h^ are 
in C^°(R d ), and by letting k — > oo we get 



h(x) > E x [h{X Tu )\ - C l6 E x [ Tu ] and h(x) < E x [h{X Tu )\ + C 16 E x [ Tu ]. 
Now, we have by (|4.21j) and (|4.27[) . for every A > 1 and x £ D(X~ 1; r, A _1 f), 

V(8 D (x)) = h(x) 

> E x [h (^r D(A _ 1?A _ lf) ) ;^r D(A _ 1?A _ 1?) € D^X-'r) \ D{\- 1 ?, X- 1 ?)] - C W E X [r D(A - 1F , A - 1?) 

> ^(A-^l + A 2 )- 1 ^^^ (x Td(a _ 1fa _ 1f) £ i?(f,A- 1 f)\^(A- 1 f,A- 1 f)) -C 16 E X [^(x-ip.A-if) 



(4.27) 



(4.28) 



We also have from (|4~27ll 

V(fe(a;)) = /i(z) < E K \h (■X'r D(x _ lftA _ 1|0 )] + Ci 6 E k [t d(a -i F , a -i F) ] 

< V(.R)P* (^r D(A _ 1? , iX _ lf) e r>) + C 16 e,[t d(a -i F:A -i F) ]. 

By (|2. 19(1 and the monotonicity of j, for every A > 4 and x £ D(X~ l r, A _1 r), 

F * (Xx*-^-^ e £>) > Px (*r D(x _ lr . x _ lF) e D(r, A- J f) \ ^(A-V, A-V)) 

(*t d(x _ 1p , x _ 1p , e ^(3A-V, A-V) \ D(2A- 1 f, X- 1 ?)) 

j(\ x s -y\)dyds 

D(3\- 1 r,\- 1 r)\D(2\- 1 r.\- 1 r) 



(4.29) 



> 
> 



T D(A- 1 ?.A- 1 P1 



27 



> / dyb'(10A- 1 r)E !C [r I)(A -i ?)A - 1?) ] 

\J D(3X- 1 r,X- 1 r)\D(2X- 1 r,\- 1 r) ) 

> c 1 (A- 1 f) d j(10A- 1 f)E ;c [r D(A -if,A-iF)] ■ 
Now, applying Theorem 12. Ill we get for x E Z?(A _1 r, A _1 r) 

Px (x TD{x _ 1?x _ lf) &D)> c 3 ^((10A- 1 f)- 2 )A a l^ [r D(A - 1F , A -i ?) ] . (4.30) 

Thus from ([£28 j) -([£30 ]) . for every x E D(X~ 1 r, \~ l ?) 

V(5 D (x)) > (^^(A-^l + A 2 )- 1 / 2 f)£((10A-V)- 2 )A Q - C 16 ) E x [t d{x -i pa -i 9} ] (4.31) 

and 

V(5 D (x)) < c 3 (l + (^((lOA-V)- 2 ))- 1 A-«) P e (x Td(x _ 1?x _ 1?) e £>) . (4.32) 

Using first (|2.12p and then Potter's Theorem ([31 Theorem 1.5.6 (1)]), we see that there exists a large 
Ao > 4 such that for every A > Ao 

U(A- 1 (l + A 2 )- 1 / 2 f)^((10A- 1 f)- 2 )A" 

>c 4 f Q / 2 (l + A 2 )- Q / 4 A Q / 2 (£((A- 1 (l + A 2 )- 1 / 2 f)" 2 )) _1/ %((10A- 1 f)- 2 )>2C 16 /c 2 (4.33) 

and 

(€((10A- 1 f)- a ))" 1 A- Q <c B . (4.34) 

Combining (|4.3ip ~ (|4.34p . we have proved the lemma with R 5 := Aq r. □ 

It is clear that every C 1,1 open set is K-fat, i.e., for any C 11 open set with C 11 characteristics (R, A), 
there exists a constant K E (0, 1/2], which depends only on (R, A), such that for each Q E dD and r E (0, R), 
D n B(Q, r) contains a ball B(A r (Q), nr) of radius nr. In the rest of this paper, whenever we deal with C 1,1 
open sets, the constants A, R and n will have the meaning described above. 

Recall that g is defined in (|3.11|) . 

Theorem 4.6 Suppose that D is a bounded C 1 ' 1 open sei m R d with C 1 ' 1 characteristics (R, A). TTien i/iere 
exists Cig = C\g(R, A, a, diam(D)) > smc/i i/iai 

Cfg 1 (V(fe(a:)) A 1) < j(x) < C 19 (V(M*)) A 1) , a: € £>, (4.35) 

or equivalently there exists C20 = C2o{R, A, a, diam(D)) > smc/j i/ioi 

, (WJL A A < , , fto 1 (Wg A A , e D . 



Proof. Since d = 1 case is simpler, we give the proof for d > 2 only. Recall that -R3 is the constant in (|3.6p 
and £i = R3K/24. Since g(x) = Go(x,Zo) A C13 and g(x) = Go(x,Zo) for 5d(x) < 6ei, it suffices to show 
that there exist r* E (0, 6ei) and ci > 1 such that 

c^ / fo ^° /a < G^,,*) < Cl /foW)^ , M.)<^. (4-37) 

v^((M*))- 2 ) " v^foWH 
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Let r* := (-R5/4) A (ei/(4>/l + (1 + A) 2 )) and suppose that Sd(x) < r* . Choose x Q G dD satisfying 
5d{x) — \x — Xq\. We fix the C 11 function ip and the coordinate system CS — CS Xo in the definition of C 11 
open set so that x = (0,x d ) with < x d < r* , 

B(x ,R) DD = {y = (y, y d ) e B(0, R) in CS:y d > i/j(y)}. 

Letxi := (0, r*/2) and D* :=D(r*,r*) = {y £ D : r* > y d ~ij(y) > 0, |j/| < r*}. Since S(xi, c 2 r*) cfl, 
for small C2 > 0, by Theorem 12.91 Theorem 12. 151 and the fact that D is bounded, 

& D (x,z ) < c 3 G D (xi,zo) 



< c 3 G(xi,z )- 



t (X TDr GB(« ,ei/4)) 
P x (X TDt eB(«b,ei/4)) 



^(^(W) efl ^ £ i/*)) 
- C4P, (X TDt g B(z ,£i/4)) < cgE, [r D J 

where in the last inequality we used (J2.19I) and the fact that dist(.D*,.B(2o,£i/4)) > 8d{zq) — £i/4 — 
-y/l + (1 + A) 2 r* > ei (see P~2l]i ). On the other hand, by Theorem [2~T5l Lemma 13731 and the fact that 
D is bounded, 

, V x (X Tn e D) . , 

G D (x,zo) > c s Gn(xi,z ) _ x \™* -4- > c 6 P x (X TOt e D) . 

Xxi i^T D , € L>) 

By applying (l4~22|) - (|433)) . we have proved ((4735]) . The inequalities (|4736)l follow from (|2~l2"|) . □ 

Now we give the proof of Theorem ll.il which is the main result of this paper. 

Proof of Theorem ll.il By [51 Theorem 1.5.3], the local boundedness and strict positivity off, there exists 
Ci > 1 such that for every < s < t < 6 diam(D) =: QM 

s a Cl t a 



Thus, if s, t, u < QM then 

_ x f s a t a u a \ (sViVuf , / s a t a u a . 

Cl ^^2) v I(Fi) v ^I^y J - ^((svtvu)- 2 ) - Cl V^F 2 ! V ^(F 5 ) V ^(tF 5 )/ ' ! ~'" kS ' 

Also note that by (13~T )) - (|3~9")) . for every ^4 € B(x,y) 



1 min Kfl5 < s < fl5 £(s 2 ) < l((£ D (A)) 2 ) < ^^ max Kfl5 < s < fis l(s 2 ) 

2 max £1 < s < A f £(s" 2 ) ~~ l{{8 D {x)\J 8 D (y)\J \x - y\)- 2 ) ~ min £1 < s < A f ^(s -2 ) 



and 



Q A j^J S D (A) < 6n(x) V 8 D (y) V\x-y\< 2k' 1 (J?- V l) S D (A). (4.40) 

By Theorem 11.21 and Theorem 14.61 we have that 

, (M*)Mj/)) q/2 *((M^))- 2 ) 1 



(Sr,(x)8r,(lA) a / 2 £((Sr,(A))- 2 ) 1 

< c 2 



^£((8 D (x))- 2 )i((5 D (y))-*)(8 D (A))<* t(\x - y\~ 2 )\x - y\ d -« ' 
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Since P~55]l - ([OP]) and the identity 



a b ab ab . „ , 

- A - A -T- = 7 : rir, for all a, b, c > 

6 a c 2 (aV6Vc) 2 



imply that 



" 3 (fefe)) a/2 Vl&WP) (M*)) q/ V*((Mj/))- 2 ) V*((MaO)- 2 )'((Mi/))- a )N-y|' 

< (fo(o ; )^(y))«/ 2 £((fe(^))- 2 ) 



vlfoM'ltfoinffofi))" 



(Ms)) a/ V*((M2/))- 2 ) (fo(^))" /2 ^((^(x))- 2 ) < (6 D (x)6 D (y)) a ' 2 l(\x-y\- 2 ) 

SC3 ; ^^^^=^^^^= A ; ^^^^=^^^^ A — " " — — 



'(MJ/)) Q/ V*((M*))- 2 ) (M*)) a/ V*((Mj/))- 2 ) Vt((SD(x))-z)£((6 D (y))-2)\x-y\*' 
to prove the theorem, we only need to show 

(M*)) Q/2 v^((M2/))- 2 ) ' {5 D {y)Y' 2 ^WM^) (5 D {x)5 D {y)Y' 2 l{\x y\~ 2 ) 



(S D (y)r/Wl((S D (x))- 2 ) (5 D (x))«/^£((S D (y)y 2 ) ^£{(5 D {x))~ 2 )£{{5 D {y))- 2 )\x - y\« 
(6 D (x)8 D (y)r/ 2 £(\x-y\- 2 ) \ 



< 1A 



^£((6 D (x))- 2 )£((S D (y)y 2 )\x ~ y\<* 



<cr (S D (x)r/ 2 ^/£((S D (y))- 2 ) a (^(2/)) a/2 y/^((^(.T))^) A {6 D (x)6 D (y)r'H(\x-v\- 2 



(fc(i/)) a/2 vf^F) (M*)) q/ V*((Mj/))- 2 ) ^((fe(x))- 2 K((<5 I) (y))- 2 )|x-yh y /- 

(4.41) 

Since the first inequality is clear, we will proceed to the second inequality. 

By symmetry, we only need to consider the cases 6d(v) < Sd(x)/3 and Sd(x)/3 < 8o(y) < 38d(x), and, 
using the fact £ is slowly varying and [2j Theorem 1.5.3], the case Sd(x)/3 < 8o{y) < 35d{x) is clear. 

Now we assume 8n{y) < Sd{x)/3. Since \x — y\ > Sd(x) — Sn(y) > 2Sd{x)/3 in this case, using [3J 
Theorem 1.5.3], the continuity and strict positivity off, 



5 D (x)6 D (y)£(\x - y 



_,.\-2 ] 2/ a 
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(*((M*))- 2 )) 1/a (*((Mi/))- 2 )) 1/o b - v\ 

< 5 D (x)6 D (y) ( (£((5 D (x))- 2 )) 1/a \ 2 Mv) (*((M*))- 2 )) 1/a 



'(<((fc(*))- 8 )) 1/a (#D(l/))- a )) 1/a l M*) / (^((^(2/))- 2 )) 1/Q feW 

and 

(M2/)) Q/ V*((M*))- 2 ) . (M*)) Q/ V*((Mi/))- 2 ) 

< c 7 - 



(M*)) q/ V*((Mj/))- 2 ) " (Mj/)) q/ V*((Mz)) 

Thus 

/ (fe(x)fe(i/)) a /^(| a: -v|- a ) ^ 



<c 8 



y/£((5 D (x))-i)£((6 D (y))-2)\x - y\<* J 
\5 D {y)r /2 ^£{{5 D {x))- 2 ) A (5 D (x)S D (y)r/ 2 £(\x - y\- 2 ) 



(Sn(x)r/^£((6 D (y))-*) ^£{{8 D {x))-*)£{{8 D {y))-*)\x-y\ 



< cg (6 D (x)r/^£((5 D (y))-2) a (fo(y))"/ 2 v Ag((M^P 2 ) a (MaOMy))"^!* - 2T 2 ) 



JiDl!/))^ 2 ^!^^))- 2 ) (M*)) Q/ V*((Mtf))- 2 ) Vt((SD(x))- 2 )£((5 D (y))- 2 )\x - y\« J 
We have obtained the second inequality in (|4.41[) . D 
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Remark 4.7 Similarly as in Theorem \4-6\ by use of (|2.12[) . the inequalities (jl.3|) imply the alternative 
forms given in (|1.4[) and (|1.6p . 



Now we give the proof of Theorem 11.31 which is a consequence of Theorems 11.11 and 12.151 

Proof of Theorem 11.31 Using the interior ball condition of D, the following holds: For every Q £ 3D 
and r < R there is a ball B — B(zq,r) of radius r such that B C D and dB n <9-D = {Q}. In addition, it 
follows from [29l Lemma 2.2] that, for each Q £ dD, we can choose a constant c 2 = C2(d, A) £ (0, 1/8] and a 
bounded C > open set [/q with uniform characteristics (i?*, A*) depending only on (i?, A) and d such that 
B(Q, c 2 i?) nflcf/QC B(Q, i?) n D and 

My) = Su a (y) for every y £ B{Q, c 2 R) n D. (4.42) 

Assume that r € (0,C2-R], Q S 5-D an d w is nonnegative function in R d harmonic in D n B(Q,r) — 
Uq n -B(Q, r) with respect to X and vanishes continuously on D c n -B(Q, r). Let zq := Zq . By [9j Lemma 
4.2] and its proof, we see that u and x — > Gu Q (x, zq) are regular harmonic in Uq n -B(Q, 2r/3) with respect 
to A. Since the C 1,1 characteristics of Uq depend only on (R, A) and d, by the boundary Harnack principle 
fTheorem I2.15|) . there exist n = ri(a,£, R, A) g (0,1/4] and c% = cs(a,£,R,A) > such that for any 
r £ (0,ri] we have 

^ Gg Q fo*Q) for every x, y ei?(g,r/2) HO. 

w(y) G[/ Q (2/,z Q ) 

Now applying Theorem 11.11 to Gjj q (x, zq) and Gu Q (y,ZQ), then using (|4.42p . we conclude that for 
r£(0,(c 2 RAr 1 )} 

u( X ) s%{xwt((6u Q iv))-*) g^ygggggj . _, n ...__ 

— r-r < c 4 — 77^ = C4 — 7^ for every x, y £ B(Q, r/2) fl L> 

for some C4 = 04(0, £, i?, A) > 0. The form (|1.7|) given in the statement of the theorem is equivalent to the 
one in the display above for r £ (0, (c 2 i? A n)]. Now the case r € {{C2R A ri), (i? A l)/4] follows from the 
case r € (0, (c 2 R A n)] and Theorem f2. 141 □ 
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